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^ . ABSTRACT 

X 

. We present a new mechanism for creating the observed cosmic matter- 

antimatter asymmetry which satisfies all three Sakharov conditions from 
one common thread, gravitational waves. We generate lepton number 
through the gravitational anomaly in the lepton number current. The 
source term comes from elliptically polarized gravity waves that are pro- 
duced during inflation if the inflaton field contains a CP-odd component. 
In simple inflationary scenarios, the generated matter asymmetry is very 
small. We describe some special conditions in which our mechanism can 
give a matter asymmetry of realistic size. 
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1 Introduction 



As far into the universe as we can see, there is an excess of matter over anti-matter. 
The recent determination of the cosmological parameters from the cosmic microwave 
background by the WMAP experiment gives the baryon density of the universe as |T] 



This is a small number, but at the same time it is large enough to be a puzzle for 
models of particle physics. A baryon excess this large cannot be produced in the 
early universe within the Standard Model of particle physics Pj. Models that can 
explain the baryon excess typically involve exotic nonstandard physics, CP violating 
couplings in the Higgs or supersymmetry sectors or in the couplings of the heavy 
neutral leptons associated with neutrino mass [Sj. In this paper, we introduce a new 
mechanism for the creation of the matter-antimatter asymmetry, one associated with 
gravitational fluctuations created during cosmological inflation. 

The conditions for generating a matter-antimatter asymmetry were stated by 
Sakharov almost forty years ago [4j. First, baryon number should be violated. Sec- 
ond, CP should be violated. Third, these symmetry violations should be relevant 
at a time when the universe is out of thermal equilibrium. Since the 1980's, it has 
been realized that the standard weak interactions contain processes, mediated by 
sphaelerons, which interconvert baryons and leptons and are thermally activated at 
temperatures greater than 1 TeV. Thus, we can also create the baryon asymmetry 
by creating net lepton number at high temperature through out-of-equilibrium and 
CP-asymmetric processes |5l6j . Scenarios of this type are known as leptogenesis. 

The out-of-equilibrium conditions can be created at a phase transition or through 
late decay of massive particles. The most attractive choice for a phase transition is 
that associated with electroweak symmetry breaking. However, that phase transition 
is probably not sufficiently strongly first-order. This is known to be an obstacle 
to baryogenesis in the Standard Model, and most of the allowed parameter space 
in minimal supersymmetry is also already excluded. Particle decay asymmetries are 
loop-suppressed and therefore require relatively large CP-violating phases. Such large 
phases are strongly constrained in supersymmetry [2j though they still could appear 
in the neutrino Yukawa couplings that are used in the Fukugita-Yanagida scenario 
for leptogenesis [6^. In any event, there is good reason to seek more effective sources 
of CP-violating out-of-equilibrium physics. 
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2 Outline of the mechanism 



We would like to assemble a theory of matter-antimatter asymmetry out of the 
following ingredients. First, as is well-known ;Sj, the lepton number current, and 
also the total fermion number current, has a gravitational anomaly in the Standard 
Model. Explicitly, 

d,Jt = jI^^RR (2) 

where 

= + , RR = ^e''^^' RafSp.Ry ■ (3) 

The anomaly requires an imbalance of left- and right-handed leptons, so we are ig- 
noring right-handed neutrinos. In general, (0) will be correct in an effective theory 
valid below a scale fi. A simple guess for is that it is at the right-handed neutrino 
scale, of the order of 10^^ GeV, but we would like to keep in mind the possibility of 
higher values of /i. 

Next, we claim, a contribution to RR of definite sign can be generated by gravi- 
tational fiuctuations produced during infiation if the infiaton field contains a CP-odd 
component. This can be naturally achieved if the infiaton is a complex modulus field 
such as one finds in supergravity or superstring models. Such fields can have the very 
flat potentials required for inflation. The imaginary part (p of this fleld (which we 
henceforth call an 'axion') can couple to gravity through an interaction 

AL = F{(f))RR , (4) 

where F is odd in (p, as a result of the Green-Schwarz mechanism [9J. Lue, Wang, and 
Kamionkowski (LWK) have studied the effects of such an interaction in generating 
observable parity- violation in the cosmic microwave background PU]- (Some other 
discussion of the effects of the RR term in gravity can be found in JI]). A simple 
form for F{(j)) with the correct scaling is 

lOTT^Mpi 

where the Mpi in the denominator is approximately the string scale. In principle, we 
could imagine substituting a lower mass scale 3^ in the denominator of (0), with the 
constraint that our effective fleld theory is valid only for yU < 5". 

We would like to apply the interaction Q to the dynamics of metric fluctuations 
during inflation. When the axion fleld has a slowly-rolling nonzero classical value, 
the coupling can lead to quantum fluctuations of the gravitational fleld that, 
treated to second order, generate a nonzero right-hand side for (0) • The scalar metric 
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perturbations generated during inflation have a size that gives density fluctations with 
6p/ p ~ 10^^. So it is interesting to speculate that a correct theory might reahze the 
relation that we see numerically, 



(6) 

Whether we succeed in realizing this intuition, we leave for the reader to judge. 



We believe that our analysis is also interesting because the Sakharov conditions 
are satisfied in this scenario in an unusual way. Lepton number is of course violated 
through (j2)). CP violation and out-of-equilibrium result from the nonzero classical 
value of the axion field. Before infiation, the complex modulus field varies from point 
to point in both modulus and phase. Infiation blows up a small region in this field to 
a size much greater than that of the visible universe. In this region, the modulus field 
is approximately constant and has a randomly chosen, fixed phase. This value then 
rolls slowly toward the minimum of its potential. In this process, we have out-of- 
equilibrium dynamics and, if the phase is nonzero, a CP asymmetry. We claim that 
no explicit CP violation is needed in the equations of motion. The CP-odd field 
could have zero expectation value today and need have no relation to the CP violation 
observed in particle physics. 



3 Gravity wave evolution 

Now we would like to quantitatively estimate the lepton number produced in 
infiation ^2]. To begin, we must compute the production of gravitational waves 
in infiation under the infiuence of the coupling (^. The general form of metric 
perturbations about an FRW universe can be parametrized as 

ds^ = -(1 + 2ip)dt'^ + Widtdx' + a^{t) [((1 + 2^)5,^ + hij) dx'dx^] (7) 

where cp, ip, Wi and hij respectively parametrize the scalar, vector, and tensor fiuctu- 
ations of the metric. It is straightforward to show that the scalar and vector pertu- 
bations do not contribute to RR, and so we ignore these fiuctuations in the following 
discussion. We can also fix a gauge so that the tensor fiuctuation is parametrized by 
the two physical transverse traceless elements of hij. For gravity waves moving in the 
z direction, we write 

ds"^ = -dt^ + a^(t) [(1 - h+)dx'^ + (1 + h+)dy^ + 2h^dxdy + dz^] , (8) 

where a{t) = e^^ during infiation and /;,+ , hy^ are functions of t, z. To see the CP 
violation more explicitly, it is convenient to use a helicity basis 

hL = ^{h+ - ih^) , hn = -^{h+ + ih^) . (9) 
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Here Hl and hji are complex conjugate scalar fields. To be very explicit, the negative 
frequency part of hi is the conjugate of the positive frequency part of hji, and both 
are built from wavefunctions for left-handed gravitons. 

The contribution of tensor perturbations to RR, up to second order in hi and hji, 

is 



~ 4i 

RR = — TT 



\dz'^ ^dtdz ^ dz^ ^ dtdz ^) \(9t^ ^ dtdz ^ dt"^ ^ dtdz ^ 



Id 2 fd , 92 g 
+ 29^^ U "^9^ ^"5^ 



(10) 



If hi^ and have the same dispersion relation, this expression vanishes. Thus, for 
RR to be nonzero, we need a 'cosmological birefringence' during inflation. Such an 
effect is induced by the addition of (jH) to the gravitational equations pUj . 

Specifically, by adding (0]) to the Einstein action, inserting (jH)), and varying with 
respect to the metric fluctuations, we find the equations of motion 

Uhi^ = -2i-h' , Dhji = +2i-h' , (11) 
a a 

where 

e = -^iF"<p' + 2HF'<P) , (12) 

dots denote time dervatives, and primes denote differentiation of F with respect to 0. 
Note that (jlj) with a constant is a total divergence that cannot affect the equations 
of motion; thus, all terms in involve derivatives of 0. We have dropped terms with 
third-order derivatives of hi and hji and terms with 0. In fact, it is also permissible to 
ignore the F" term in ()12|). since in slow-roll inflation, <^ Mp\H and each derivative 
on F brings a dimensionful factor of order 1/Mpi. These evolution equations should 
be compared to those for evolution in flat space given by LWK pi3|. The new term 
proportional to H(j) leads to a substantial enhancement in the size of 0. With this 
simplification, and the approximate form (0), 

= v^K(iJ/Mpi)V27r2 , (13) 

where e = ^(0)^/(ifMpi)^ is the slow-roll parameter of inflation |12j . 

Let us now focus on the evolution of and, more specifically, on its positive 
frequency component. It is convenient to introduce conformal time 
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(Note that conformal time t] runs in the opposite direction from t.) The evolution 
equation for becomes 



hL-2-—hL--f^hL = -2te-—hL 15 

If we ignore G for the moment and let hi ~ e*'^^, this becomes the equation of a 
spherical Bessel function: 

^hL-2-^hL + k''hL={) (16) 

drj'^ T] drj 

for which the positive frequency solution is 

hl{k,7]) = e+*^(''+")(l -ifcr/) . (17) 



We now look for solutions to (|T5|) with /il ~ e**^^. To do this, let 

Hl = e'^' ■ (-ifcr7)e^%(r7) 
Then girf) satisfies the equation 
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drf^ 



k\i _ e^) - ^ - ^ 



. (19) 



This is the equation of a Schrodinger particle with £ = 1 in a Coulomb potential. 
When B = 0, the Coulomb term vanishes and we find the spherical Bessel function 
ffTTj) . For hi, the small Coulomb term is repulsive; for hji, with the opposite sign of 
the G term, the Coulomb potential is attractive. This leads to attentuation oih^ and 
amplification of in the early universe. This is just the cosmological birefringence 
described by LWK pUj. 

It will turn out that the generation of the matter asymmetry is dominated by 
modes at short distances (sub-horizon modes) and at early times. This corresponds 
to the limit kr] ^ 1. In this region, we can ignore the potential terms in and 
take the solution to be approximately a plane wave. More explicitly, 

g{7]) = exp[tk{l - 62)1/2^(1 + a{7]))] , (20) 

where a (77) ~ log rj/rj. 
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4 Green's function 



We would like to apply the results of the previous section to compute the expec- 
tation value of RR in the inflationary space-time. Our expression will be dominated 
by the quantum part of the gravity-wave evolution. For this regime, we can calculate 
the expectation value by contracting and in RR using an appropriate Green's 
function. 

Define 

-^^e^H^-^'^Gk{v,v') ■ (21) 

For k parallel to z, the Fourier component Gk satisfies (I15|) with a delta- function 
source 

ar]"^ rj arj arjj Mp^ 

For B = 0, the solution of this equation is 



Gk{v,v)=^^-T7^S{v-v') . (22) 



UHy2k'M^,)ht{k,v)h],{-k,r^'), r^<ri 
^oyV,V) ]^^Hy2k^M^,)hlik,v)hU-k,v') v' < V ^ ' 

where h~[ is the complex conjugate of (fTTj) . and h^, are the corresponding solutions 
of the hn equation. For 9 = 0, these solutions are the same as for /i^, but the structure 
of will be preserved when we go to the case 9 7^ 0. The leading effect of 9 is to 
introduce the exponential dependence from (fTHj) . 

The prefactor is modified in order 9^, and the wavefunctions acquire additional cor- 
rections that are subleading for krj ^ 1. Neither of these effects will be important 
for our result. 



The Green's function (j24j) can now be used to contract hi and to evaluate the 
quantum expectation value of RR. The result is 



im = - / 7^ T^TJTwikvr ■ k'Q + 0(9^) (25) 



16 r d^k H"^ 
T J (27r)3 Wm^^ 

where we pick up only the leading behavior for kr] ^ 1. 



We note again that our expression for (RR) is nonzero because of the effect of 
inflation in producing a CP asymmetry out of equilibrium. The original quantum 
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state for the inflaton might have had nonzero amphtude for a range of values of 
and might even have been CP-invariant. However, inflation collapses the wavefunction 
onto a particular value of (j) that is caught up in the local expansion of the universe. 
This value gives us a classical background that is CP-asymmetric. 



5 Lepton asymmetry 

Now we are ready to evalute the lepton density that arises through ((21). Inserting 
fl25|) into (j2I) and integrating over the time period of inflation, we find for the net 
lepton number density 



The integral over k runs over all of momentum space, up to the scale jj, at which our 
effective Lagrangian description breaks down. The dominant effect comes not from 
the usual modes outside the horizon at the end of inflation (super-horizon modes), 
k/ H < 1, but rather from very short distances compared to these scales. The integral 
over 1] is dominated at large values of rj, early times. The integral represents a 
compromise between two effects of inflation, first, to blow up distances and thus 
carry us to smaller physical momenta and, second, to dilute the generated lepton 
number through expansion. It is now clear that the dominant contribution to the 
right-hand side comes from kr] ^ 1, as we had anticipated. Performing the integrals, 
we find 



We might interpret this result physically in the following way. The factor (if/Mpi)^ 
is the usual magnitude of the gravity wave power spectrum. The factor B gives the 
magnitude of effective CP violation. The factor is the inverse horizon size at 
inflation; this gives the density n appropriate units. Finally, the factor (fi/H)^ gives 
the enhancement over one's first guess due to our use of strongly quantum, short 
distance fluctuations to generate RR, rather than the super-horizon modes which 
effectively behave classically. 

To understand the significance of this estimate, we should compare it to the en- 
tropy density of the universe just after reheating, assuming that the energy of the 
inflationary phase has been converted to the heat of a gas of massless particles. To 
estimate this, assume very naively that reheating is instantaneous. Then reheat- 
ing converts an energy density p = SH^Mp^ to radiation with p = n'^g^.T'^/SO and 




(26) 




(27) 
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s = 27r^5f*T^/45, where (yf* is the efffective number of massless degrees of freedom. 
This gives s = 22g^'^{HM^if/^ . With this value 
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Assuming that there have been no large increases in the entropy of the universe 
since the end of reheating, ()28|1 can be compared directly to the present value of 
n/s inferred from (0). For this one should note that the ratio of the present baryon 
number to the lepton number originally generated in leptogenesis is approximately 
n^/n = 4/11 0; then (P) implies n/s = 2.4 x 10"^°. 

Inserting the estimate for B given in ()13p and setting ~ 100, we find 

The final result has n/s (x corresponding to n oc H\ I gnoring for the moment 

the dependence on /x, the various factors accumulate to a very strong dependence on 
{H/Mp\). Unfortunately, the ratio (H/Mpi) is limited in simple slow-roll inflation from 
the relation Sp/p ~ {H/M-pij/^/t ~ 10^^. The WMAP results give a more precise 
version of this bound for the case of single-field inflation: H/Mpi < 1 x 10~^ fTS]. To 
make an estimate, put /i ~ if ~ 2 x 10^^ GeV; this yields 

n/s < IQ-^^ ■ v^N , (30) 

a very small result for any reasonable choice of JnT. 



6 Improvement of the lepton asymmetry 

The estimate ()3()|1 is rather unsatisfactory, but in several places we have used 
naive or pessimistic values of parameters. We should now ask by what factor (jHn|) 
could reasonably be enhanced. Looking back to (PHjl . the prefactor and the factor 
(if/Mpi)''/^ are intrinsic to our mechanism. The remaining factors, however, are 
negotiable. 

In the estimates above, we assumed that the mass scale from the kinetic term for 
the modulus is the Planck or string scale. If this mass scale is set at a lower mass 3", 
9 can be larger, scaling as 

e ~ V2^y^{H/Mpi)\Mpi/3')/27i^ . (31) 
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In principle, 'N can be a large dimensionless number, though values greater than 100 
are probably difficult to accomodate with string theory. In the estimate (IHUj). 6 is of 
the order of 10~^. It is straightforward to make 6 larger by three to four orders of 
magnitude. 

Another possible source of enhancement is already apparent in ()28|1 . We can 
obtain a much larger value of n/s by raising the value of fi at which our effective 
Lagrangian is valid. Taking 3^ fi and setting g^, ~ 100, K ~ 100, and e ~ 10~^ we 
obtain 



H is bounded above by the WMAP data, but there is also a lower bound on H from 
the fact that in a leptogenesis model sphaelerons must be active after inflation. This 
requires a reheating temperature T > 1 TeV, corresponding, in a model with fast 
reheating, to H > 10~^ eV. For H in this range, the observed value of n/s is found 
for 



The upper end of this range is somewhat higher than the conventional scale of super- 
symmetric grand unification. The lower end is within the conventional range for the 
masses of right-handed neutrinos. For if ~ 1 TeV, fi ~ 10^^ GeV. 

In summary, we have presented a new mechanism for the production of the cosmic 
baryon asymmetry. This mechanism relies on the axial vector anomaly to violate 
fermion number and on the initial state of inflation to both for CP violation and 
for out-of-equilibrium dynamics. These are very minimal ingredients that might be 
found in a wide variety of models of physics at short distances. The simple version of 
this model that we have analyzed here generates an acceptable values for the matter 
asymmetry only if we can make use of gravitational fluctuations at scales of the order 
of the usual grand unification scale. It is interesting to ask whether the conditions 
we have found can be embedded in a grand unification model in a natural way. 
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